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Abstract. If an axion of mass between about 10 -3 eV and 10 eV exists, axion emission 
would have significantly affected the cooling of the nascent neutron star associated with 
SN 1987A. For an axion of mass less than about 10 -2 eV axions produced deep inside 
the neutron star simply stream out; in a previous paper we have addressed this case. Re- 
markably, for an axion of mass greater than about 10 -2 eV axions would, like neutrinos, 
have a mean-free path that is smaller than the size of a neutron star, and thus would 
become “trapped” and radiated from an “axion sphere.” In this paper we treat the “trap- 
ping regime” by using numerical models of the initial cooling of a hot neutron star that 
incorporate a diffusion approximation for axion-energy transport. We compute the axion 
opacity due to inverse nucleon- nucleon, axion bremsstrahlung, and then use our numerical 
models to calculate the integrated axion luminosity, the temperature of the axion sphere, 
and the effect of axion emission on the neutrino bursts detected by the Kamiokande II 
(KII) and Irvine-Michigan-Brookhaven (IMB) water-Cherenkov detectors. The larger the 
axion mass, the stronger the trapping and the smaller the axion luminosity. We confirm 
and refine the earlier estimate of the axion mass above which trapping is so strong that 
axion emission does not significantly affect the neutrino burst: Based upon the neutrino- 
burst duration — the most sensitive “barometer” of axion cooling — we conclude that for an 
axion mass of greater than about 3 eV axion emission would not have had a significant 
effect on the neutrino bursts detected by KII and IMB. The present work, together with 
our previous work, strongly suggests that an axion with mass in the interval 10 -3 eV to 3 
eV is excluded by the observation of neutrinos from SN 1987A. 
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I. Introduction 

Peccei-Quinn (PQ) symmetry may be the simplest and most compelling extension of 
the standard SU(Z)c 0 SU(2)l 0 £f(l)y model. PQ symmetry cures the single blemish 
on QCD: the strong-CP problem, and predicts the existence — but not the mass — of a 
new pseudoscalar particle: the axion. 1 A priori the mass of the axion could be anywhere 
between about 10 -12 eV and 1 MeV, corresponding to PQ symmetry breaking scales 
between about 10 19 GeV and 100 GeV. (The axion mass and PQ symmetry breaking scale 
are related by m a /eV ~ 6 x 10 6 GeV/(/ a /iV); the axion coupling to ordinary matter is 
proportional to m a — or equivalently, (J a /N)~ l .) A host of astrophysical and cosmological 
arguments — and a few laboratory searches — have left open but two “windows” for the 
axion mass: 10 -6 eV to 10 -3 eV and about 1 eV to 5 eV (hadronic axion only); see 
Refs. 2 for an up-to-date review of the “axion window.” One of the most powerful and 
important constraints to the axion mass is based upon the early cooling of the neutron 
star associated with SN 1987A. Axion emission can accelerate the cooling of the nascent 
neutron star and thereby shorten the neutrino burst. In particular, it has been argued that 
the neutrino bursts detected by the Kamiokande II (KII) and Irvine- Michigan- Brookhaven 
(IMB) water- Cherenkov detectors would have been significantly shorter than the bursts 
actually observed if an axion in the mass interval of 10 -3 eV to 2 eV existed. 3 (Many 
authors have studied the possible effect of axions on the cooling of SN 1987A; Ref. 4 
contains a semi-complete bibliography.) 

At the temperatures and densities relevant to the hot, newly bom neutron star, the 
dominant process for axion emission (and absorption) is nucleon- nucleon, axion brems- 
strahlung (and inverse axion bremsstrahlung): N + N <-► N + N + a. Axion emission from 
the nascent neutron star can be divided into two qualitatively different regimes: “freely 
streaming,” for m a < 0.01 eV; and “trapping,” for m a > 0.01 eV. In the freely stream- 
ing regime the axion-mean-free path for absorption is large compared to the size of the 
neutron star, and axions, once emitted, simply “freely stream” into the vacuum of space. 
In the trapping regime, axions interact sufficiently strongly that their mean-free path for 
absorption is small compared to the size of the neutron star; in this case, like neutrinos, 
they are said to be “trapped” and are effectively emitted from an axion sphere. (The 
axion sphere is the surface beyond which the probability for an axion to be absorbed is 
exp(— 2/3) ~ 0.5.) 

Neglecting the “back reaction” of axion emission on the cooling of the neutron star, 
axion emission in the freely streaming regime is simply proportional to the axion-nucleon 
coupling squared which is proportional to the axion mass squared. In the trapping regime 
things are more complicated; in the simplest treatment, the axion luminosity is propor- 
tional to the fourth power of the temperature of the axion sphere. Based upon a simple 
analytic model 3 (which this work shows to be quite good) it has been argued that the tem- 
perature of the axion sphere varies as m a 4 ^ X1 , so that the axion luminosity in the trapping 
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regime should vary as m a 16/,n . Very roughly then, one expects that as a function of axion 
mass, the axion luminosity should increase as m 2 for m a <C 10 -2 eV, and should decrease 
as ma 16/,X1 for m a >• 10 -2 eV (see Fig. 1). From this simple picture, one sees that here 
should be two “critical” masses for axion emission from SN 1987A: one below which axion 
emission is acceptable because the axion interacts so weakly; and one above which axion 
emission is again acceptable because the axion interacts so “strongly.” 

The freely streaming regime is relatively simple to treat: A heat sink of magnitude 
equal to the local axion-emission rate is incorporated into numerical models of nascent 
neutron star cooling. In previous work we did just that. 5 Based upon the duration of the 
neutrino bursts that would have been observed in the KII and IMB detectors we concluded 
that the “lower mass boundary” is about 10 -3 eV . (Several other studies are in agreement 
with our conclusion. 4 ) The trapping regime is more difficult to address because in principle 
one has to treat axion-energy transport in much the same way as one does neutrino-energy 
transport (or radiative transport in an ordinary star). Based upon a simple analytic 
model the “upper mass boundary” was estimated to be about 2 eV. 3 The existence of the 
previously mentioned axion window around a few eV depends crucially upon the upper 
mass boundary: Were it 5 eV rather than 2 eV the window would be closed. Moreover, 
one experiment to search for axions in this mass range is currently being carried out, 6 
and another has been proposed. 7 The first involves searching for the photon-line radiation 
produced by the decays of relic (cosmological) axions; 6 and the second involves detecting 
axions emitted by the sun by axion-photon conversion induced by a strong magnetic field. 7 
For this reason, and the general importance of the SN 1987A bound to the axion mass, we 
are addressing axion transport and emission in the trapping regime. 

To preview our results, the window doesn’t “close up.” Based upon the present work 
we conclude that the upper boundary mass is about 3 eV, rather close to the original 
estimate of about 2 eV. The present work together with our previous work 5 strongly 
suggests that the durations of the neutrino bursts detected by KII and IMB exclude an 
axion with mass in the interval 10 -3 eV to 3 eV. We axe quick to remind the reader that 
both mass boundaries depend upon the precise form of the axion-nucleon coupling, as well 
as the neutron star models and the exact burst-duration exclusion criterion. The mass 
boundaries are therefore “fuzzy” by about a factor of two or so. 

The outline of the rest of the paper is as follows. In the next Section we will calculate 
the crucial physics input to the problem: the axion opacity (under ordinary circumstances, 
this would be an oxymoron). In Section III we will derive the equations that govern axion 
transport and the diffusion scheme that we employ. Section IV is devoted to a discussion 
of the results of our numerical simulations of axion- cooled neutron stars, and in Section V 
we summarize and add some concluding remarks. 
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II. Axion Opacity 

As we have discussed above, for an axion mass of greater than about 10~ 2 eV, it 
is expected that the axion-mean-free path for absorption (at densities and temperatures 
typical of a newly born neutron star) is less than the radius of a neutron star. 3 In this mass 
regime axions do not simply stream out and one has to calculate the axion luminosity in 
much the same way one does the photon luminosity in an ordinary star or the neutrino 
luminosity in a newly bom, hot neutron star. To do so one needs to calculate the axion 
opacity as a function of density, p, temperature, T, and the axion energy, E a . The axion 
opacity, he, at energy E a is related to the axion-mean-free path by 

(kep) -1 = K(E a ,p,T). (1) 

In the present circumstance, unlike photon transport in an ordinary star or neutrino 
transport in a hot neutron star, only absorption is important. This is because each 
axion line in a Feynman diagram introduces a dimensionless coupling factor of order 
m iv/(/o/A) ~ 10 -7 (m a /eV), and so processes involving more than one axion are sup- 
pressed relative to those involving a single axion by a factor of at least 10 -14 (m a /eV) 2 . 
By far the dominant axion-absorption process is inverse nucleon- nucleon, axion brems- 
strahlung (a + A r + A — ► iV + iV; A is a neutron or proton). 

There are two equivalent methods for computing A a . The first, more familiar to a 
physicist, relies upon the Boltzmann equation. The mean-free path for absorption is related 
to the attenuation of the phase space density of a stream of axions moving in the z direction: 

1 = fa 1 ~qZ~ = 2E~ j ^ 1< ^ 2( ® 3C ^ 4 (^ 7r ) 4 ^ 4 (P a + Pi + P'2 ~ P 3 — P 4 ) 

x5|M| 2 / 1 / 2 (1-/ 3 )(1-/4), (2) 

where pi, p 2 , £> 3 , P 4 and p a are the nucleon and axion four momenta, the subscripts 1, 2 (3, 
4) refer to the incoming (outgoing) nucleons, dll, = d 3 p t /(2ir)*2E t is the Lorentz-invariant 
phase-space- volume element, /,• are the nucleon-phase-space distribution functions, 5 is 
the symmetry factor (a factor of 1/2 for identical nucleons in the initial or final states), 
and \M\ 2 is the matrix-element squared (summed over initial and final spins). Throughout 
this Section we shall set h = kg — c = 1. For reference, we remind the reader that the 
axion-nucleon interaction follows from the Lagrangian density 

Ant h (garx/2m N )(n'f li 'y 5 n)d >l a + (flfa P /2m JV )(p7 M 7 5 p)5 /i a, 

where the axion-nucleon couplings g an — c n m^/(f a /N) and g ap = c p m^/(/ a /A), and c p 
and c n are numerical constants of order unity. For more about the axion and its couplings 
to ordinary matter see Refs. 8. For the derivation of the matrix element squared and the 
details of carrying out the phase-space integrations see Ref. 9. 
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The second method for computing A a , more familiar to an astrophysicist, relies upon 
Kirchhoff’s law (also known as detailed balance or time-reversal invariance) for calculating 
the opacity, 

E dp a (T)/dE' ( } 

where A a is related to Kg by Eq. (1), j e is the axion emission rate (at energy E) per gram 
of material per second per axion energy interval, and 


1 El 

dE. 2tt 2 e E 'l T — 1 ’ 


is the differential axion energy density for a thermal distribution of axions. (Expression 
(3) is probably even more familiar when written for photons: j„ = An k vcl B v {T), where 
B U (T) = 2hv 3 / (e hv t kT — l)c 2 is the Planck function, dp^/dv — AnB v jc, and u = E/h = 
E/2n is the frequency; also note that because axions are spinless particles, the Planck 
functions for photons and axions differ by a factor of 2.) The total axion-volume- emission 
rate, i a (used in previous work on axions and SN 1987A), is related to Je by 

f ° ° 

£a = p 3EdE a , (5a) 

JO 



*S\M?hh(l - /,)( 1 - /,)( 1 + /.). (56) 

To make the calculation of A a tractable we will make some approximations. First, 
we will assume that |Af| 2 is approximately constant; as discussed in Refs. 10 this is a 
reasonable approximation at the temperatures and densities of interest (we will have more 
to say about this below). Assuming that g a = g an = g ap , 5|Af| 2 is given by 9 

S|A4| 2 = ^P^-(3-/?) (6a) 

3 my 

for n + n + a— >n + n and p + p + a— +p + p, and 

5|A<| 2 = ^^(7-2/3) (66) 

for n + p + a—*n + p, where a n = (/m^/m f ) 2 ~ 56 is the nucleon-pion coupling factor 
and /? is a parameter that depends upon the degree of nucleon degeneracy. For completely 
degenerate nucleon matter /? — * 0, and for non-degenerate matter j3 — » 1.0845; see the 
Appendix of Ref. 9 for further details. Next, we assume that the nucleons can be treated 
as being non- degenerate. Deep in the core this is a marginal approximation; 9 however 
further out, near the axion sphere (T ~ lOMeV, p ~ 10 12 g cm -3 ) where all the action 
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is, this is a good approximation. (Ishizuka and Yoshimura 10 have recently computed A a 
in the degenerate limit.) We also assume that the nucleons are non-relativistic, which is 
a very good approximation throughout the star. (In Ref. 11 the fully-relativistic matrix 
element and phase-space integrations are compared to the non-relativistic matrix element 
and phase-space integrations.) Finally, in the most important region, that near the axion 
sphere, the densities are such that many-body effects, e.g., reduction in the effective nu- 
cleon mass and variation of the pion-nucleon and axion-nucleon couplings, should not be 
significant (see Ref. 11). In sum, the ambient conditions near the axion sphere are such 
that the various approximations we make are well justified. 

With these approximations it follows that: 


/. = exp(t/j - u,-), Ui= p 2 j2m N T, y,- = (m - m N )/T, 
n i = — ^==(m N T) 3/2 e yt , 

7rv27r 

where rij is the number density of nucleon species i (= n or p ) and /i, is the chemical 
potential of species i. With these approximations the axion-volume-emission rate for the 
process N\ + N 2 —* + N 4 + a is given by 9 

= 4 ^5^5 m N 2 T 6 ' 5 e Mvi + ya). 

It is now straightforward to evaluate analytically expression (2) for A" 1 : 


,_i _ S\M\ 2 nin 2 

“ ~ 2 5 tt 3 / 2 m 5 /2Ti/2 ( 


,EJT 


i)l. 


du-uj 2 (u- - E a /2T) 1/2 e ~ 2u - . (7) 


E a /2T 


(The integrals in expression (2) for A" 1 are evaluated in an analogous manner, and using 
the same notation, as those for i a are in Ref. 9.) The final integral factor which is a 
function of E a /T can be expressed in terms of the modified Bessel function A'i(x): 12 


e E./T(r\ /•“ du_uL /2 (u- - E a /2T) lf2 e~ 2u ~ 

\EJ Jb./1T 

1 1 T / P \ 

= -exp(£ a /2T)A' 1 (£ 0 /2T)~ i — (l + ^) , (8) 

where the second expression is a useful empirical fit which has an accuracy of about 
10%. Note that the axion-mean-free path is relatively insensitive to the axion energy and 
temperature, and varies roughly as 


A a oc T 1 / 2 


EJT 

( 1 + Ejry / 2 ’ 


( 9 ) 
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this is in agreement with the original estimate made in Ref. 3. Note too that the inverse 
of the axion-mean-free path is proportional to the target-number density squared, rather 
than the target-number density; this of course is because the absorption process has a 
three-body initial state with two target nucleons. 

Because the axion is a boson, the presence of ambient axions will lead to stimulated 
emission of axions, cf. the factor of (1 + f a ) in Eq. (5). Owing to this fact, the “net” 
absorption (= true absorption less stimulated emission) is less than the true absorption cal- 
culated above, and a “reduced” absorption opacity is often defined. Assuming an ambient 
thermal distribution of axions, the reduced absorption opacity is 


k* e = k e (l - e E ‘/ T 'j . 


( 10 ) 


The quantity k* e describes the net axion absorption as a flux of axions passes through 
matter. As one can readily see the reduced absorption opacity and the absorption opacity 
do not differ by a large factor since the typical axion energy E a ~ 3T. It is also useful to 
define the Rosseland-mean opacity 




&Pa(T) 

dEdT 



&Pa(T} 

dEdT 


dE, 


( 11 ) 


which weights k e near the peak of the energy flux (E a = AT). Using the energy dependence 
of the axion- absorptive opacity computed above, cf. Eq. (7), we find numerically that 



( 12 ) 


In order to compute the total axion opacity one must consider all three absorption 
processes (n + n + a — > n + n, p + p + a — * p + p, and n + p + a — * n + p); this is 
accomplished by adding the corresponding expressions for A^ 1 from each process: 


A a x (total) = A a 1 (nn) + A a \pp) + A a : (np); (13) 


opacities, like the resistances of resistors in series, add. Finally, taking g an = g ap = g a = 
2 m N/(fa/N) ~ 7.6 x 10 -8 (m a /eV) and setting (3 = 1.0, we have evaluated A~ x (total) 
numerically: 




(m a \ 2 

( t r /2 i 

r p \ 

1 u 

VMeVy ' 

^10 14 gcm -3 J 


x (1 + $,X n X p )e Ell / 2T Ki(E a /2T); 

~ (2.4 x lO’cm)- (^) 2 ( ( ,/ A 

VeV/ \MeV/ \10 14 gcm~ 3 y 


2 


2 


(14a) 
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( 144 ) 


X (1 + SX.X r ) (£) (l + fr)'^; 

where X n and X p are the neutron and proton mass fractions, and in Eq. (14b) we have 
used our empirical expression for exp(E a /2T)Ki(E a /2T). We have taken g an = g a p = 
/(fa/N) for consistency with previous work; 5 the axion-nucleon couplings depend 
upon the type of axion — DFSZ or hadronic — the quark distribution functions that are 
assumed (because of the 75 couplings, it is the fraction of the nucleon spin carried by quarks 
that is relevant— and still uncertain). As we will remind the reader in our concluding 
remarks, our results depend upon the assumed values of the couplings — and can be re- 
scaled for different assumed values of the axion-nucleon couplings. From Eqs. (14) we can 
compute the axion opacity at energy E a 


Kg == 2.1 x 10 


-18 


cm 2 \ /m a \ 2 f T \ 1/2 


)W ( 


Kg ~ 4.2 x 10 


-18 


g ) VeV 
x (1 + 8X n X p )e F ‘/ 2T AV(£ a /2T); 

2 \ \ 2 / T \ -1 / 2 


10 14 gcm 


-3 


cm 


g 


VeW \MeV J 


10 14 gcm -3 
1/2 


xa + «jr,)g)(i + ^) . 

Finally, we can compute the Rosseland-mean opacity; using Eq. (12) we find: 


M* “ 195 * 10-18 (?) OY (i Sv)~ I/2 (1 + 8W (15) 


(15a) 


(156) 


Of the approximations made in calculating the axion opacity — non-degenerate and 
non-relativistic nucleons, and constant matrix element squared — the latter is least well 
justified. Because of the various pion-propagator factors that enter in the matrix element 
squared, |A4| 2 is not constant. The dependence upon the nucleon-momentum transfer 
enters in the form of the following pion-propagator factors: 


l*'| 4 /(|£| 2 + ml) 2 , |/|V(|/i 2 + m 2 ) 2 , |fc| 2 |/i 2 /(|fc| 2 + mlm 2 + m 2 ); 


, 2\2 


see Ref. 9. Here m n = 135 MeV is the pion mass, and k = p 2 — and l = p 2 — P3 
are the four-momenta transfers in the two types of Feynman diagrams (for details, see 
Ref. 9). The three-momenta exchanged are |£| 2 ~ |/| 2 ~ 3 m^T; at high temperatures, 
T m 2 /2m n — 6 MeV, the pion-propagator factors become momentum independent and 
equal to unity. Deep inside the core, the temperatures are sufficiently high that the pion- 
propagator factors can be ignored; further out — say near the axion sphere where T ~ 10 
MeV — the validity of ignoring the pion-propagator factors is less justified. 
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To be more quantitative about the effect of the pion propagator we have computed the 
axion-emission rate with a pion-propagator factor included, j P g , by multiplying the con- 
stant matrix element squared by |£| 4 /(|£| 2 + m 2 ) 2 , and comparing it to the rate computed 
without the pion-propagator factor — the canonical assumption. (We note that this proce- 
dure is not precisely correct, as the pion-propagator factor that occurs in the interference 
terms in |A4| 2 involves both |fc| and |/| — see Ref. 9; however, this procedure should give 
one a pretty good idea of the effect of the pion propagator.) 

The ratio of the pion-propagator corrected rate to the uncorrected rate is given by 

Tf( T \ _ r e_I K x + <*) 1/2 z 1/2 ( 1 + 4e 2 /y+y_) - elm (y+/y-)] dx 
} j~t-*{x + otyiWidx 

where a = E a /T, e = ^ and y± = [(x + a) 1 / 2 ± x 1 / 2 ] 2 + 2e. The axion-emission 

reduction factor R is shown as a function of temperature in Fig. 2 for E a /T = 1, 4, 7, and 
10. For E a /T = 4 (a characteristic value for a thermal distribution and for the Rosseland 
mean) and T > 10 MeV, R is greater than about 1/2. Since k# oc jg, R is a measure of the 
reduction in both axion emission and axion opacity due to the effect of pion propagators 
in the matrix element squared. 

It is straightforward to show that R has the following limiting behaviors: R — * 1 — 0(e) 
for e — + 0 and R — > e~ 2 for e 1. Motivated by this we have used the following expression 
to approximate R: 

R(a, T ) = 7— T — T-y. 

1 + a(a)e + b(a)e 2 

For a = E a /T = 4, a = 0.814 and b = 0.054 give the correct limiting behaviours and a fit 
that is accurate to better than 7% for all values of T. (For reference, a = 1.027, b = 0.0673 
for a = 3 and a = 2.22, b = 0.107 for a = 1.) Although the effect of including the pion 
propagators is small (see end of Section IV), we have used the above fit to R for a = 4 to 
correct both i a and k e for use in the numerical models in Section IV. 

III. Axion-Energy Transport 

To properly treat the effects of axions upon the cooling of the nascent neutron star 
associated with SN1987A in the trapping regime ( m a > 10 -2 eV) one must employ the 
full apparatus of radiative- transfer theory. This is a formidable task, and in light of all the 
uncertainties involving in setting a limit to the axion mass — neutron star equation of state, 
the initial state of the hot neutron star, and the imprecision of our “exclusion criterion” 
based upon the length of the neutrino burst, we have opted to use an approximation 
scheme — diffusion — to describe the transport of axion energy out of the newly born neutron 
star. 

Axion transport in SN 1987A is very similar to photon (radiative) transport in an 
ordinary star, and so we will adopt the language and machinery that has been developed 
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for that problem. (For the most part we follow the notation and conventions used in 
Refs. 13, except we use the axion energy E instead of the axion frequency v.) The primary 
quantities of interest axe the various moments of the specific intensity, Ie- The specific 
intensity describes the flow of energy carried by axions ( dE a ) in a particular direction (n) 
through an area (dA) into a solid angle (dQ) per energy interval (dE = hdv) per time (dt), 

dEg 

E cos QdAdSldEdt' 


where 0 is the angle between dA and n. The specific intensity Ie is related to the axion 
phase space density f a by Ie = E 3 f a /h 3 c 2 . 

The equation of radiative transport (which follows directly from the Boltzmann equa- 
tion) governs the evolution of Ie and is given by 


1 dls 

c dt 


+ n • VIje = 


~P*e(Ie ~ Be), 


(17) 


where local thermodynamic equilibrium (LTE) is assumed for the nucleon species, Be 
is the previously discussed Planck function (for axions), and only axion-absorption and 
emission processes have been included. Taking the zeroth and first angular moments of 
this equation, specializing to spherical symmetry, and integrating over axion energy we 
obtain 

^^ +p “^T 1+ 5^l (4,rr2 - F - ) = - c / - E(dp./dE-4*B E /c)dE , (18) 


F 

® a 


C dpg . 

3 p{k)r dr F ' 


(19) 


Here p a = J(dp a /dE)dE is the axion-energy density, p a is the axion pressure, F a is the 
radial axion-energy flux, and we have used the Eddington approximation, p a ~ 3p a , which 
is valid for any nearly isotropic radiation field (i.e., everywhere, except very near the axion 
sphere). 

The general-relativistic version of Eq. (18) was solved by an operator-split method, 
along with the neutrino-energy transport equations, and the equation of hydrostatic equi- 
librium. A major simplifying assumption was to take the axions to be in LTE with the 
matter interior to the axion sphere. In the strongly trapped regime this is a very good ap- 
proximation. In the transition regime (total “axion-optical depth” r = dr / \ a between 

1 and 100), corresponding to axion masses between 2 x 10 -2 eV and about 0.2 eV, this is 
a marginal approximation as the axion sphere is broad and ill defined. As we discuss in 
the next Section we were unable to treat the transition region reliably for other reasons 
too. However, we are confident that the the effect of axion cooling in the transition is so 
severe that such an axion mass is clearly excluded. 

Finally, we mention that we previously tried a less robust approximation to treat axion 
transport: leakage. In the leakage approximation the effect of axion emission is modeled 
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by a local heat sink of magnitude £„/( 1 + kr 2 ) where k is a numerical constant of order 
unity and r is the total axion optical depth. While leakage schemes have been successfully 
used in other instances, e.g., neutrino transport, 14 they were not up to the difficult task 
of treating axion-energy transport. This is in part because the leakage approximation is 
well suited to the case where the mean-free path and density axe constant; in the present 
circumstance the density varies significantly near the axion sphere (and A a oc p -2 ). The 
leakage scheme systematically underestimated the effect of axion cooling. 

TV. Results of Numerical Models 

For the purpose of this investigation of the trapped regime, we have focused on a single 
proto neutron-stax model, model B from our previous work on the freely streaming regime 
(Ref. 5, hereafter referred to as BTB). In BTB, we investigated a variety of proto-neutron 
star models, with different equations of state and different masses, and found substantially 
the same results for all of the models. We feel confident that we can restrict the present 
studies of the trapping regime to a single model, model B, the model that best reproduces 
the neutrino observations of SN 1987A. (Model B has a stiff equation of state and the 
proto-neutron star mass starts at 1.3 M© and increases by accretion to about 1.5 M©.) 

For the most part we will follow the approach of BTB. To briefly remind the reader 
of the strategy of the previous work, we first computed the neutrino flux from a numer- 
ical simulation of the cooling of the nascent neutron star that included freely streaming 
axion emission. From this flux and the response characteristics of the KII and IMB wa- 
ter Cherenkov detectors we computed the number of i/ e -capture events expected for each 
detector ( N ) and the time required for the expected number of events to reach 90% of 
its asymptotic value (Ai 90 %), again for each detector. In addition, we computed the total 
energy carried off by axions ( E a ) and by neutrinos (E„). 

Both the energy carried off by neutrinos and the number of capture events were only 
mildly sensitive to the effects of axion emission; we found that the most sensitive indicator 
of axion emission was A t 90 %. As the assumed axion mass was increased to about 10 -3 
eV, the neutrino burst duration dropped precipitously to less than about 1 sec for both 
detectors, in contrast to the observed burst durations of about 12 sec (KII) 17 and about 
6 sec (IMB); 18 see Fig. 3. On this basis, in BTB we concluded that the KII and IMB 
data excluded an axion more massive than about 10 -3 eV (at least in the freely streaming 
regime). 

To investigate the trapping regime we used the same general-relativistic code employed 
in BTB, which is discussed in detail in Refs. 15 and 16. Numerical simulations of the initial 
cooling phase (first 20 sec after collapse) of a hot neutron star that incorporate the axion 
diffusion were carried out for axion masses from 0.5 eV to 10 eV. For each model we 
computed the total energy radiated in axions E a , the total energy radiated in neutrinos 
E„, and for each detector, the expected number of events N and burst duration A f 90 %. 
In addition, at times of 0.01 sec, 0.1 sec, 1.0 sec, 5.0 sec, and 20 sec we computed the 


11 



temperature and luminosity of the axion sphere (the bulk of the axion luminosity comes 
from the axion sphere), red shifted to r = oo. Our results are given in Tables I, II, and III 
and illustrated in Figs. 1 and 3-6. In addition, both for completeness and comparison, in 
Figs. 1, 3, and 4 we have shown the results of our previous work for the freely streaming 
regime. 

For model B the total axion-optical depth becomes of order unity for an axion mass 
of about 0.02 eV, in good agreement with the estimate previously given in Ref. 3. If for 
the moment we neglect the effect of axion cooling on the structure of the hot neutron star, 
then the total axion-optical depth scales simply as r ~ (m„/0.02eV) 2 . From this we see 
that axions are not strongly trapped (say r > 100) unless m a > 0.2 eV. In this awkward 
transition regime, m a ~ 0.02 eV to 0.2 eV, the axion opacity is significant, and axions are 
not freely streaming; however, they are not so strongly trapped that there is a well defined 
axion sphere. Neither treatment is suitable for the transition regime. Moreover, in the 
transition regime the axion luminosity is near maximum — recall, the simplest expectation 
is that the axion luminosity peaks for r ~ 1 — so that the effect of axion cooling is very 
great. 

In attempting to treat the transition regime we discovered another difficulty: During 
the initial collapse from the white dwarf state {p ~ 5 x 10 9 gem -3 and T ~ 0.5 MeV) and 
just before the proto neutron star state is reached, the total axion optical depth will reach 
unity and the axion luminosity will be very large (10 55 ergs sec -1 in one case we studied). 
So large in fact that the initial model we assume for the nascent neutron star is not self 
consistent as the collapsing core will certainly be affected significantly by axion emission 
during the collapse. In order to self consistently consider an axion mass of between 0.02 eV 
and 0.2 eV one would have to take into account the effect of axion emission on the collapse 
phase. For these reasons we decided that the transition region was beyond the scope of 
the present investigation, and only considered axion masses greater than 0.5 eV, for which 
r ~> 1000 and the effect of axion emission on the collapse and initial state of the neutron 
star should not be very large. We would argue strongly that the effect of axion cooling 
for axion masses between 0.01 eV and 0.5 eV is so severe that an axion of this mass is 
most certainly precluded. For such masses there is every reason to expect that the energy 
carried away by axions is greater than that for a mass of 0.5 eV, and therefore that the 
burst duration should be shorter. We should emphasize that because of the dependence of 
the axion-mean-free path on density and temperature, for any axion mass between about 
0.03 eV and 30 eV axions will, at the onset of collapse, be freely streaming and then, only 
late in the collapse, will they become trapped. As they become trapped their luminosity 
will be achieve its maximum value. Thus we expect that axions of mass greater than about 
0.02 eV will affect the collapse and the initial state of the proto neutron star. 

Because axion emission during the initial collapse phase can be significant — and we 
do not take it into account — we decided that it was important to check the consistency of 
our treatment. Our chief concern is the fact the the initial state we assume for the proto 
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neutron star is not quite “relaxed.” That is, because energy transport by axions was not 
included during the collapse phase, the temperature gradient and lepton number gradient 
are not quite what they should be; had we included the effect of axion cooling during the 
collapse phase, the initial proto neutron star would be relaxed. To gauge the importance 
of the “mismatch” in initial conditions, in Table I we show the energy radiated in axions 
before and after the proto neutron state relaxes. For an axion mass of 0.5 eV these energies 
sue almost equal; for a mass of 1 eV the energy radiated before relaxation is only about 
1/3 of the total energy radiated; and and for a mass of about 3 eV it is only about 25%. 
We should emphasize that the mismatch in the initial state only leads to an uncertainty 
in the energy radiated while the neutron star is relaxing to equilibrium. Owing to this 
fact we estimate an uncertainty in the predicted number of neutrino events of about ±1/2 
event which, for our purposes is not significant. 

The behavior of our best barometer for axion emission — the burst duration At 90 % — is 
shown in Fig. 3 for axion masses from 10 -4 eV to 10 eV, spanning both the freely streaming 
and trapping regimes. One can clearly see the effect of axion emission on the neutrino burst: 
As the axion mass is increased to about 10 -3 eV the burst duration drops precipitously 
due to the effect of axion emission; at a mass of about 10 -2 eV, where trapping begins, 
the burst duration is less than 1 sec for both KII and IMB. For masses between 0.01 eV 
and 0.5 eV we strongly suspect that the burst duration remains shorter than about 1 sec, 
but as discussed above we have been unable to reliably treat this regime. At a mass of 0.5 
eV the IMB burst duration is still less than about 0.5 sec and the KII burst duration is 
less than 2 sec. As the axion mass is increased beyond 0.5 eV trapping reduces the effect 
of axion emission, and A t 90 % increases. For an axion mass of 2 eV the predicted burst 
duration for KII is about 5 sec, about half of the predicted duration in the absence of axion 
emission. For an axion mass of about 4 eV the IMB burst duration has increased to 2 sec, 
about half of the predicted burst duration in the absence of axion emission. On this basis 
we would set the upper boundary of the excluded mass interval to be about 3 eV. 

The reader is reminded that the physical origin of the precipitous drop in burst dura- 
tion, first found in BTB, can be traced to the fact that there are two distinct phases of 
neutrino emission. The first phase, lasting of order 1 sec, is powered by the heat in the 
outer mantle and residual accretion; the second phase, lasting of order many seconds, is 
powered by the outward diffusion of the heat trapped in the inner core. The first phase 
is rapid because the timescales for neutrino diffusion out of the low density outer mantle 
and for residual accretion are both short (< 1 sec). The second phase lasts much longer 
because the timescale for diffusion of neutrinos from the inner core is long, of order many 
seconds. Axion emission tends to deplete the heat trapped in the inner core that powers 
the second phase of neutrino emission by providing another means of transporting heat 
out of the inner core. By so doing, axion emission can drastically shorten the duration of 
the neutrino burst. 

As in the freely streaming regime, the number of capture events (see Fig. 5) and the 
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energy carried off by neutrinos (see Fig. 1) are much less sensitive indicators of axion 
emission: At a mass of about 10~ 2 eV the expected number of events falls by less than a 
factor of 2, and, at most, axions carry away as much energy as the neutrinos do. This is 
simple to understand: Axion emission does not directly suppress neutrino emission; rather, 
axions tap the same source of energy as do neutrinos, and thus, axion cooling serves mainly 
to shorten the cooling time. 

There is one new interesting twist in the trapping regime: The number of events 
expected for the IMB detector rises above its value for no axion cooling at an axion mass 
of 0.5 eV, before relaxing to its no axion cooling value for large axion mass. This odd 
phenomenon has a simple explanation: For axion masses around 0.5 eV, axions are very 
efficient in transporting energy from deep in the core out to the neutrino sphere, thereby 
heating the neutrino sphere. Because of the high energy threshold of the IMB detector, 
the number of events expected at IMB is very sensitive to the temperature of the neutrino 
sphere. 

The total energy carried away by axions is shown in Fig. 1 and Table I. In the regimes 
where axions axe a minor heat sink, axion masses much smaller than 10“ 2 eV or much 
greater than 1 eV, one can, for purposes of understanding how E a scales, ignore the back 
reaction of axions on the cooling of the neutron star. Doing so, in the low mass regime 
one expects the energy carried off by axions should vary as m 2 since the axion luminosity 
is proportional to i a which varies as m 2 . In the high mass regime, the situation is more 
complicated because of axion trapping. However, one expects the axion luminosity to vary 
as the temperature of the axion sphere to the fourth power, and in Ref. 3 it was estimated 
that the temperature of the axion sphere should vary as this implies that the 

energy carried off by axions should vary as rria 16/fn in the large mass limit. Both of these 
scalings are roughly consistent with our numerical results for the temperature of the axion 
sphere and the energy carried away by axions. 

It may be of some interest to know the average energy of the axions emitted by a 
nascent neutron star in the trapping regime, e.g., if one envisions nascent neutron stars 
as an intense source of axions that might be detected by other means. 19 In Table II and 
Fig. 5 we have shown the temperature of the axion sphere at times of 0.01 sec, 0.1 sec, 
1.0 sec, and 20 sec. (The axion sphere is the surface beyond which the number of axion- 
mean-free paths, or “optical depth,” equals 2/3; i.e., J^p(/c)^dr = f£?dr/\ a = 2/3.) 
The temperature has been red shifted to r = co: T 0 0 — T a /(1 + z) where (1 + z) is the 
gravitational red shift from the axion sphere to infinity. (The total energies carried off by 
axions and neutrinos discussed earlier were red shifted in the same way.) The average axion 
energy is related to the temperature of the axion sphere by (E a ) = 2.7T a /(l + z). One can 
see that the temperature of the axion sphere drops roughly as a power law, Too oc f -1 / 4 , 
from which it follows that the energy radiated in axions per logarithmic interval of time 
is roughly constant. In Fig. 6 we display the axion luminosity from the axion sphere at 
the same times; very roughly, C a decreases as f -1 , as one would expect since C a oc T 4 and 
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T a cx i" 1 / 4 . 

Next, we briefly comment upon the inclusion of the pion-propagator correction factor 
R(a,T). In general, its effect upon our results was small (typically 10% to 20%) and 
would not have significantly affected any of our conclusions (had we not included it). As 
expected, the effect of including this correction was most significant for large axion masses, 
for which the temperature of the axion sphere is the smallest (recall that T a cx m a 4 ' /11 ). 
Based upon the very small effect of including this correction, we feel confident that not 
including the pion-propagator in our previous work 5 was a well justified approximation. 
In particular, we re-ran model B for freely streaming axions of mass 10~ 2 eV, including 
the pion-propagator correction, and the results changed insignificantly. 

Finally, we wish to illustrate the richness involved in adding the effect of axions to the 
initial cooling of a hot neutron star. For comparison, we have run some initial cooling 
models without axion cooling, but with the addition of extra light (mass 10 MeV) 
neutrino species. In Figs. 7 and 8 we show the predicted number of events and burst 
duration A t 90 % for the KII and IMB detectors as a function of the number of neutrino 
flavors, N„ = 3 to 11. (The additional neutrino species, like the ^ and r neutrinos were 
assumed to couple only through the neutral current interaction.) As can be seen from 
these Figures, the effect of additional neutrino species has a much less dramatic effect: 
The number of events and burst duration decrease slowly with N„. The reason is simple. 
Since all neutrino species couple with weak interaction strength, their effect on transport 
of heat from the inner core is minimal. (Very roughly, the effect of additional neutrino 
species on neutrino energy transport is to increase the effective neutrino coupling strength 
by a factor of (iVy/3) 1 / 2 .) The additional species serve mainly to dilute the energy carried 
off by electron antineutrinos, because the energy released from the gravitational collapse 
of the core must be shared among more degrees of freedom. In contrast, as the axion 
mass is varied the coupling strength varies (as m a ). By varying m a from 10 -4 eV to 
10 eV, one explores a range of qualitatively different regimes, from freely streaming to 
strongly trapped, and the effect on the cooling of the newly bom neutron star is much 
more dramatic. While one can exclude an axion in the mass range of 10~ 3 eV to about 
3 eV based upon the KII and IMB observations, based upon the same observations one 
would be hard pressed to exclude as many as 5 additional neutrino species. Of course, 
the existence of more than 3 light neutrino species is now firmly excluded by the precise 
measurements of the width of the Z° boson made by the SLC experiment at SLAC and 
the ALEPH, DELPHI, L3, and OPAL experiments at LEP. 19 

V. Concluding Remarks 

The existence of an axion of mass in the range 10 -3 eV to about 10 eV would have 
had a significant effect upon the cooling of the nascent neutron star associated with SN 
1987A (and upon other newly born neutron stars). For such a mass, axions would carry 
away a significant fraction of the energy and would significantly accelerate the cooling 
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process. We have now investigated in detail the effect of axion cooling both in the freely 
streaming regime ( m a < 10 -2 eV) and in the trapping regime ( m a > 10 -2 eV). Based 
upon the duration of the expected neutrino bursts calculated in our axion-cooled models, 
an axion mass in the interval 10~ 3 eV to 3 eV can be excluded. The upper mass boundary 
is very close to the original estimate made in Ref. 3. Therefore, the axion mass window 
around a few eV remains the same, from about 3 eV to 5 eV (hadronic axions only). The 
entire window was explored by a search for the decays of relic cosmological axions that was 
carried out at Kitt Peak this year. 6 If that search is unsuccessful, it will close the multi-eV 
axion-mass window. The multi-eV window is also accessible through an experiment that 
has been proposed to search for axions emitted by the Sun. 7 It is also possible that this 
mass region could be explored if a supernova explosion occurred within our own galaxy, 
e.g., by more closely examining the neutrino signal (provided that many more events are 
detected) or by other means such as gamma-ray observations. 20 

As we emphasized earlier in this paper and in BTB, our results, which are expressed in 
terms of the axion mass, actually depend upon the values of the axion-nucleon couplings; 
for definiteness we have assumed that c n = c p = | where g an = c n m_\r / (f a /N) and 
g ap = c p m^/(f a /N). The dimensionless axion-nucleon couplings c n and c p depend upon 
the PQ charges of the quark species and the quark-distribution functions; the couplings 
are discussed in some detail in Refs. 21. Both mass boundaries for the excluded mass 
region scale with the inverse of the axion-nucleon couplings; that is, doubling c n and c p 
would decrease both the upper and lower mass boundaries by a factor of two. (Of course, 
the re-scaling of the boundaries of the excluded region is more difficult if c„ and c p do not 
change in the same way; however, one could probably still estimate the change.) 

Finally, we should mention the uncertainties inherent in our axion mass constraint. 
To begin with, there are the uncertainties associated with our numerical cooling models — 
equation of state, neutron star mass, amount of residual accretion, the diffusion approxi- 
mation used for axion transport, and our exclusion criterion for the duration of the neu- 
trino burst. While these uncertainties could amount to a factor of two or so in the mass 
boundary, additional uncertainty beyond that does not seem likely. The uncertainty in the 
axion-emission rate and opacity are a different matter. Deep in the core of the nascent 
neutron star the densities certainly reach several times that of nuclear matter; there may 
be high-density effects, nuclear many-body effects — or even an exotic form of matter at 
the core, e.g., quark matter or a pion condensate — that could significantly affect the ax- 
ion emission rate or opacity. The high-density effects have discussed in Ref. 11 and do 
not seem likely to affect either the axion luminosity or opacity significantly. In any case, 
the uncertainties associated with the high densities at the core of the neutron star would 
probably only affect the low mass boundary, as only in the freely streaming regime does 
most of the axion luminosity come from the core. In the trapping regime, most of the 
axion luminosity comes from the axion sphere, which around the upper mass boundary is 
characterized by rather modest temperatures, around 10 MeV, and densities, around 10 13 
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g cm 3 ). 

To be more specific, it was pointed out in Ref. 11 that axion emission from the core 
could be strongly suppressed if the core was quark/gluon matter rather than hot nu- 
clear matter, since the axion-quark coupling is significantly smaller than the axion-nucleon 
coupling. Very recently, Ellis and Salati 22 have considered this effect quantitatively and 
concluded that the existence of a quark/gluon matter core would significantly raise the 
lower boundary of the excluded region. Based upon the present work, we would argue 
that the lower boundary would not be raised to more than about 0.01 eV, the mass where 
trapping sets in. For an axion mass up to about 0.5 eV, the axion luminosity during the 
initial collapse phase (before the quark/gluon core can form) is prohibitively large. Fur- 
thermore, in the trapping regime, axions are in thermal equilibrium from the axion sphere 
(densities around 10 13 gem -3 ) inward to the boundary of the quark/gluon core; thus, the 
absence of a thermal bath of axions deep in the quark/gluon core should not significantly 
affect the axion luminosity, which arises primarily from axion production near the axion 
sphere. (Of course, we are simplifying matters; the presence of a quark /gluon core would 
also greatly lessen the effect of axion energy transport from the core outward — and in some 
cases inward. Firm conclusions await a detailed treatment of proto neutron star cooling 
models with a quark/gluon core.) 

Almost since its conception, it has been realized that the axion could significantly 
affect the cooling of stars of all kinds. Because the evolutionary timescales for most stars 
are measured in millions, if not billions, of years, the astrophysical arguments based upon 
stellar evolution that have been used to constrain the axion mass have necessarily been 
indirect. 2 The lone exception is SN 1987A; here the 19 neutrino events detected by KII 
and IMB provide the complete cooling history of the newly bom neutron star. Based upon 
that cooling record an axion mass in the range 10 -3 eV to 3 eV is excluded. Not only 
is this constraint the most stringent astrophysical constraint to the axion mass, but the 
directness of the argument is most pleasing. 
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FIGURE CAPTIONS 


Fig. 1. Total energy carried off from SN 1987A by axions (solid curve) and neutrinos (broken 
curve) as a function of axion mass. The results for m a < 10 -2 eV were taken from 
our previous work 5 (model B). The results for m a > 0.5 eV are from the present work. 
As explained in the text for technical reasons we were unable to consider axion mass 
between 0.03 eV and 0.5 eV; however, we do not expect any surprises in this mass 
interval. In agreement with simple arguments, the energy carried off by axions scales 
very roughly as m 2 for small axion masses and as m~ l G for large axion masses. 

Fig. 2. The pion-propagator reduction factor, R(a, T), as a function of temperature for a = 
E a jT = 1, 4, 7, and 10. The pion-propagator reduction factor is the factor by which 
the pion propagators in the matrix element for nucleon- nucleon, axion bremsstrahlung 
reduce both the axion-emission rate and opacity relative to the approximation where 
the pion-propagator factors are ignored, cf. Eq. (16). 

Fig. 3. The expected neutrino-burst duration, Af 90 %, in the KII and IMB detectors for axion- 
cooled nascent neutron star models as a function of axion mass. The results for m a < 
10 -2 eV were taken from our previous work 5 (model B). The results for m a > 0.5 eV 
are from the present work. The quantity Af 99 % is the time required for the expected 
number of neutrino events to achieve 90% of its asymptotic value. As explained in the 
text for technical reasons we were unable to consider axion mass between 0.03 eV and 
0.5 eV; however, we expect that the burst duration is very short in this mass interval. 
Based upon the expected duration of the neutrino burst, axion masses in the interval 
10 -3 eV to 3 eV are excluded by the KII and IMB data. Note that the burst duration 
in both detectors is reduced by about the same factor. 

Fig. 4. The expected number of neutrino-capture events for the KII and IMB detectors for 
our axion-cooled neutron star models as a function of axion mass. The results for 
m a < 10~ 2 eV were taken from our previous work 5 (model B). The results for m a > 0.5 
eV are from the present work. As explained in the text for technical reasons we were 
unable to consider axion mass between 0.03 eV and 0.5 eV. The intriguing increase in 
the number of events predicted for IMB for axion masses around a few 0.1 eV is due to 
the fact that axion energy transport heats the neutrino sphere, and because of its high 
threshold the IMB detector is very sensitive to the temperature of the neutrino sphere. 
As described in Section IV, the expected number of events is relatively insensitive to 
the effect of axion cooling. 

Fig. 5 The temperature of the axion sphere (red shifted to r = oo) at times of 0.01 sec, 0.1 
sec, 1.0 sec, 5.0 sec, and 20. sec. 

Fig. 6. The axion luminosity from the axion sphere (red shifted to r = oo) at times of 0.01 
sec, 0.1 sec, 1.0 sec, 5.0 sec, and 20. sec. 

Fig. 7. Predicted burst duration Atg 0 % as a function of the number of light (mass <C 10 MeV) 
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neutrino flavors. These models do not incorporate axion cooling and are shown only 
to compare the effect of additional neutrino flavors with the effect of axion cooling. 
Fig. 8. Predicted number of events as a function of the number of neutrino flavors. These 
models do not incorporate axion cooling. 
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Table I. Summary of axion-cooled nascent neutron star models for axion masses from 0.5 
eV to 10 eV. All models were run for the first 20 seconds after collapse. 


Axion 

mass 

(eV) 

Number 
of events 

KII IMB 

Energy (10 51 ergs) 
Axions Axions vv 

(before) (after) 

A*90% 

(sec) 

KII IMB 

0.5 

9.0 

3.7 

40 

53.7 

147.5 

1.7 

0.4 

1.0 

9.1 

7.4 

21 

44.9 

160.3 

3.5 

0.6 

2.0 

9.9 

7.0 

12 

28.2 

181.3 

5.0 

1.0 

3.0 

10.2 

6.5 

7 

19.8 

193.0 

6.S 

1.6 

4.0 

10.2 

6.0 

5 

15.5 

198.0 

7.5 

2.2 

5.0 

10.2 

5.8 

3 

12.4 

202.0 

S.O 

2.7 

7.0 

10.2 

5.3 

2 

8.6 

207.4 

8.7 

3.6 

10.0 

10.3 

5.1 

1 

5.9 

213.3 

9.1 

4.1 
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Table II. Temperature of the axion sphere, red shifted to r = oo and given in MeV, as a 
function of time. 


m a (eV) 

10 ms 

100 ms 

1.0 sec 

5.0 sec 

20.0 sec 

0.5 

15.1 

9.7 

5.6 

3.9 

1.7 

1.0 

13.2 

9.5 

4.6 

3.6 

2.0 

2.0 

11.3 

7.2 

4.2 

3.2 

1.95 

3.0 

10.2 

6.0 

3.S 

2.8 

1.7 ■ 

4.0 

9.6 

5.7 

3.S 

2.7 

1.7 

5.0 

9.4 

4.9 

3.5 

2.6 

1.6 

7.0 

3.9 

4.4 

3.0 

2.3 

1.4 

10.0 

7.9 

3.7 

2.6 

1.9 

1.2 
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Table III. Axion luminosity from the axion sphere, red shifted to r 
10 51 ergs sec -1 , as a function of time. 

= oo and given in 

m a { eV ) 

10 ms 

100 ms 

1.0 sec 

5.0 sec 

20.0 sec 

0.5 

490 

65 

6.6 

1.5 

0.058 

1.0 

370 

71 

33 

1.2 

0.1 

2.0 

215 

28 

2.4 

0.73 

0.095 

3.0 

1S4 

15 

1.75 

0.45 

0.06 

4.0 

157 

14 

1.67 

0.40 

0.055 

5.0 

160 

8.3 

1.29 

0.30 

0.045 

7.0 

150 

5.85 

0.72 

0.19 

0.026 

10.0 

124 

3.0 

0.46 

0.10 

0.012 
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